Diffusion of liquid domains in lipid bilayer membranes 
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We report diffusion coefficients of micron-scale liquid domains in giant unilamellar vesicles of 
phospholipids and cholesterol. The trajectory of each domain is tracked, and the mean square 
displacement grows linearly in time as expected for Brownian motion. We study domain diffusion 
as a function of composition and temperature, and measure how diffusion depends on domain size. 
We find mechanisms of domain diffusion which are consistent with membrane-dominated drag in 
viscous L Q phases [P.G. Saffman and M. Delbriick, PNAS 72, 3111 (1975)], and bulk-dominated drag 
for less viscous L a phases [B.D.Hughes et al., J. Fluid Mech. 110, 349 (1981)]. Where applicable, 
we obtain the membrane viscosity and report activation energies of diffusion. 

PACS numbers: 68.35.Fx, 68. 55. Ac, 87.16.Dg 



Diffusion of domains within cell membranes is a highly 
relevant biophysical problem. The presence of lipid do- 
mains, including rafts, can affect both short-range (intra- 
domain) and long-range (inter-domain) diffusion of mem- 
brane components PJ3l- Diffusion has been observed in 
live cell membranes [3J, |j] , although interpreting results 
from these complex systems can be challenging. 

Even in simple model systems, deciphering the diffu- 
sion of membrane inclusions is a long-standing and dif- 
ficult hydrodynamic problem [H, 0, 0] . Objects that dif- 
fuse in the membrane plane range from small peptides 
and individual lipids to large inclusions like protein ag- 
gregates and lipid domains. The first challenge is that, 
unlike in three dimensional (3D) diffusion, the size of the 
diffusing object is not the only length-scale that enters 
into the problem. For example, the membrane has finite 
thickness, a finite surface area, and often a nonzero cur- 
vature. Secondly, the membrane is composed of macro- 
molecules, which limits continuum approaches to large 
objects. Lastly, the membrane and its surroundings have 
different viscosities. In complex biological membranes, 
additional length-scales may be important, such as the 
distance between membrane proteins Q or the size of 
corrals created by the actin cytoskeleton 

In this Letter, we directly measure diffusion of liquid 
domains in giant unilamellar vesicles (GUVs) of radius 
~ 20/im as in Figure [T] These domains are micron-scale, 
circular, span the lipid bilayer, and undergo Brownian 
motion. By measuring diffusion of bilayer domains over a 
wide parameter range of more than one decade in domain 
radii and three decades in 2D membrane viscosities, we 
probe the two limiting models of Saffman-Delbruck [B| 
and Hughes et al. Q. We find a cross-over between the 
two models which would not have been predicted from 
previous monolayer results [ill ]. In the cases where our 
data are well fit by the Saffman-Delbruck equation, we 
are able to extract viscosities of lipid phases and diffusion 
activation energies. 

Domains move in a background phase with two dimen- 



sional (2D) membrane viscosity (?/'). The diffusion coef- 
ficient of a membrane inclusion was originally described 
by Saffman and Delbriick [5J: 

where r is the radius of the inclusion, 7 = 0.5772 and we 
have chosen boundary conditions appropriate for liquid 
domains in a liquid membrane to yield the factor of 1/2. 
A key parameter in the hydrodynamics of this system is 
the lengthscale Ao defined by the ratio of the membrane 
77" to the 3D bulk viscosity of water (rj w ) such that Ao = 
■q" /rj w . Eq. [I] is expected to hold for r < Ao, i.e. small 
domains and/or large membrane viscosity. Later work 
confirmed this calculation and derived a result for the 




FIG. 1: Fluorescence microscopy phase diagram of 
DOPC/DPPC/cholesterol and vesicle images at 20° C. Semi- 
quantitative dashed tie-lines cross the L a -L coexistence re- 
gion Some vesicles studied have a continuous L a (bright) 
phase (a-b) whereas others have a continuous L (dark) phase 
(c-d). One composition (e) has a continuous dark L phase 
which may contain both L and gel (S ) phase lipids [l(J. 
Vesicle compositions are shown as mol% DOPC/DPPC/Chol. 
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FIG. 2: Greyscale fluorescence images (left) are filtered and 
thresholded (middle). White regions are identified as do- 
mains. Those within a specified size range, ellipticity and 
distance from the edge are retained. Circles (right image) 
identify those domains successfully identified through 5 suc- 
cessive frames. Mean square displacement data for domains 
with radii of 1-1.5/nn are shown at right. Both vesicles have 
composition 1:2 DOPGDPPC + 30%Chol, T=10°C (top) 
and 20° C (bottom). Note the factor of ten difference in dif- 
fusion coefficients. The scale bar is 40/im. 
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It is important to notice that for Xq < r the diffusion 
coefficient is more strongly dependent on the inclusion's 
radius, but is independent of the membrane viscosity. 
This case was verified experimentally through observa- 
tions of micron-scale domains in monolayers [11] ■ Other 
theoretical work has addressed different inclusion shapes 
as well as large domains 12|. 

Spherical giant unilamellar vesicles (GUVs; 30-100/im 
diameter) are made by electroformation [l3| of a ternary 
mixture of cholesterol with phospholipids of both high 
(DPPC; di(16:0)PC) and low (DOPC; di(18:l)PC) melt- 
ing temperatures. Materials and methods have been de- 
scribed previously [Toj . The vesicle membranes are ini- 
tially uniform at high temperature, and phase separate 
into two liquid phases when vesicle suspensions are placed 
on a pre-cooled microscope stage. The less viscous L a 
phase is labeled by fluorescent dye (Texas Red-DPPE). 
The composition and viscosity of the two phases depend 
on the composition and temperature of the entire vesi- 
cle. With time, domains coalesce, allowing us to probe a 
range of domain sizes at constant temperature. 

We probe five lipid compositions in the ternary sys- 
tem of DOPC/DPPC/Chol, as in Figure Q] In previous 
microscopy and 2 H NMR measurements, we established 
that vesicles with these five compositions separate into 
a liquid-ordered (L Q ) phase rich in the saturated lipid 
DPPC and a liquid crystalline (L a ) phase rich in the 
unsaturated lipid DOPC 0, Two of the five com- 
positions contain a continuous bright L Q phase (Figure 
QJi-b), and two contain a continuous dark L phase (Fig- 
ure [2>d). One composition falls within a three phase re- 
gion (FigureQJ.) The presence of three phases is clear in 



2 H NMR experiments (manuscript in preparation), but 
is difficult to detect by microscopy. We probe the viscos- 
ity of the continuous phase by tracking domains of the 
minority phase. 

Membrane domains are identified by an image process- 
ing program written in Matlab (Figure [J) . A Gaussian 
filter is applied to images before thresholding, to identify 
features in the size range of domains. Almost no domains 
are lost by this algorithm. Domains are accepted if: a) 
the diameter falls between a minimum (2 pixel) and max- 
imum value; b) the shape is circular, such that all points 
in the domain perimeter lie within 20% of the mean do- 
main radius; and c) the center of mass lies within a circle 
defined by 0.8 of the vesicle radius. Since all domains 
are round (shape fluctuations are minimal), these crite- 
ria discriminate against occasional problems arising from 
image analysis filtering (for example, two domains very 
close to each other will not be accepted). 

A separate program tracks domain trajectories with 
logic similar to existing codes, i.e. by matching a domain 
with the nearest feature in the next image [15j . Average 
diffusion is subtracted to yield unbiased domain motion. 
Matching generates no false positives but does not have 
a perfect success rate. We therefore divide each movie 
(typically 100 frames at 0.34 s/frame) into 20 sets of 5 
frames over which most domains are tracked successfully. 
Domain size does not change over this period. The aver- 
age of vertical and horizontal mean square displacements 
(MSD) is linear with time t and fit to < x 2 >= 2D(r)t 
as expected for diffusion. Over five frames, the MSD is 
< 0.03^tm 2 which is much smaller than the particle sep- 
aration, and we see no effect of domain packing. 

Figure[3]shows diffusion coefficients as a function of do- 
main size. The data has been culled to report only sets 
in which wide ranges of domain sizes are observed for any 
fixed temperature and composition. The dashed lines in 
Figure[3]show the asymptotic 1 jr behavior given in Eq.[2] 
Within error, all data fall on or below this theoretical up- 
per bound in diffusion coefficient. Eq.[2]is independent of 
membrane viscosity and holds when membrane viscosit; 
is low, or domain radius is large i.e. when Ao < r 
For the low viscosity L a phase (e.g. Figure[3ji), it can be 
seen that the conditions of low membrane viscosity and 
large domain radius are met through most of the tem- 
perature range, because most data fall along the dashed 
line. 

All data below the dashed line in Figure [3] correspond 
to membranes with high viscosity, notably L Q phases at 
low temperatures. We have chosen to fit our data to 
the Saffman-Delbruck equation, which should hold when 
membrane viscosity is high. It is clear that this set of 
data does not have a D(r) ~ r _1 dependence, and we 
find instead reasonable fits to Eq. [T] with a single fit- 
ting parameter (?/'). Since our domain radii are limited 
to >0.5^m by our optical resolution and <10/mi by our 
vesicle diameters, we cannot prove that Eq. [T]is the only 
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FIG. 3: Diffusion coefficients vary with domain radius. Solid 
lines show fits to a logarithmic dependence on domain size 
as in the Saffman-Delbriick Eq. [1] These fits have only one 
free parameter, as discussed in the text. Symbols identify 
increasing temperatures, recorded each plot. Dashed lines 
are fits to Eq. [4] with no free parameters. Error bars report 
standard deviations for sets with >3 measurements. 



expression that could fit our data. Nevertheless, fitting 
to Eq. Q] allows us to extract 77" for high viscosity mem- 
branes. 

There have been multiple experiments designed to 
test the logarithmic form of Eq. [1] 0, [HI, llfij . and its 
range of applicability is still controversial. For exam- 
ple, recent work asserts that individual proteins diffuse 
with a stronger size dependence, as D(r) ~ r _1 
due to a break-down of the continuum approximation 



TEGS. Here we 



of the membrane for small inclusions 
explore domain radii well within the continuum limit 
(r > Single molecule ~ 0.5nm). 

Our results can be extrapolated to estimate the diffu- 
sion rate of raft domains in cell membranes. Lipid rafts 
are reported to have diameters of 10-100nm [18[. This 
length-scale falls within the regime where the Saffman- 
Delbriick equation should apply. We calculate diffusion 
coefficients for 10-100nm domains in our system to be 
between 3 x 10~ 3 and 1.5 x 10 _1 /im 2 /sec. These val- 
ues differ from those extrapolated from single molecule 
measurements using the Saffman-Delbriick equation [![. 
Given the suggestion that single molecule measurements 
do not probe the continuum limit, it may be more valid to 
estimate raft diffusion coefficients by extrapolating down 
from large domains rather than up from single molecules. 

Our analysis of the culled data set in Figure [3] shows 
that high viscosity membranes produce data that fall well 
below the dashed line in Figure [3] and that fit Eq. Q] rea- 
sonably well. We conclude that any remaining unplotted 
data that fall well below the dashed line should also fit 
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1:1 DOPC/DPPC + 30%Chol ♦ 248±2 
1:2 DOPC/DPPC + 30%Chol • 252±2 
1:4 DOPC/DPPC + 20%Chol ▼ 123±8 

FIG. 4: The diffusion coefficient Do and the 2D membrane 
viscosity rj" as a function of temperature in membranes with 
majority L and S -L phases. Solid lines are fit to Do ~ 
exp(— Ea/ksT), with E a values shown in the table. 



Eq. [TJ We use this data to yield a size-independent Do 
using: 
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Figure Q] shows a plot of the resulting Do and 2D mem- 
brane viscosity vs. temperature (T). We find that the L Q 
phase viscosities for membranes of 1:1 DOPC/DPPC + 
30%Chol and 1:2 DOPC/DPPC + 30%Chol are similar, 
suggesting that L Q viscosities are not highly composition 
dependent. In contrast, viscosities for membranes of 1:4 
DOPC/DPPC + 20%Chol are much larger, consistent 
with these membranes falling within the three phase re- 
gion in Figure [TJ 

At 22°C, we find 2D membrane viscosities for the L Q 
phase of 10~ 8 < 7/'(Ns/m) < 5 x 10~ 7 . Surface shear 
rheometry finds monolayer viscosities on the order of 
10~ 8 to 10 _6 (Ns/m) only in liquid condensed phases [lij ]. 
which is consistent with tight packing of lipids in bilayer 
L„ phases. In contrast, the 2D membrane viscosity of the 
L a phase is low, and results in D^ r _1 . In monolayers, 
the same D^ r _1 dependence is found for solid domains 
diffusing across a liquid background 11 1 . 

In the membrane literature, a 3D membrane viscos- 
ity, tizd, is defined as r]3D — Tj" '/h, where h is the bilayer 
thickness. Assuming h = 3.3nm, we find 3 < rfap(Pas) < 
150, on the order of or greater than [21( pub- 

lished values for model membranes. However, the re- 
lation between r]3D and r)" is not exact, because lipids 
anchored to the interface differ from a thin homogeneous 
layer. Indeed, the lipid headgroups often determine the 
membrane viscosity This is not always appreciated, 
and may be a source of ambiguity in discussions of the 
Saffman-Delbriick model [16J. 
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Figure [4] demonstrates that domains diffuse in mem- 
branes of high viscosity via an activated process. If the 
2D membrane viscosity (77") were independent of T, we 
would expect D (T) ~ T. Instead, we find a better fit 
for log(D (T)) ~ — T -1 , consistent with an activation 
energy E a for diffusion such that D ~ exp(— Ea/ksT). 
The data in Figure 4 follow Arrhenius behavior even 
though a gel phase emerges at low temperature for some 
mixtures. Composition of the L D phase varies only 
slightly with temperature [lfl ] . Activation energies for in- 
dividual lipids 22, H, 24 1 have been attributed to the en- 
ergy required to hop into an available free volume @, HiJ • 
Larger particles such as protein aggregates yield lower 
apparent activation energies [l7| . Fig. |4] lists activation 
energies for domains diffusing in L G phases. We find ac- 
tivation energies greater than those reported for single 
molecules in similar membranes, including DPPC/Chol 
membranes at high temperature (30-80kJ/mol) [22j, as 
well as L c lipids in phase separated DOPC/DPPC/Chol 
membranes at low temperature (~ 80kJ/mol) 

In summary we present a simple method for quanti- 
fying the movement of domains in membranes with co- 
existing liquid phases. We find that domains diffuse via 
Brownian motion, and that diffusion rates are described 
by different models under different experimental condi- 
tions. At high temperatures and in membranes with a 
continuous h a phase, membrane viscosity is low, diffu- 
sion constants are independent of membrane properties, 
and domains diffuse with a radial dependence of D ~ 1 jr. 
In membranes with a higher viscosity continuous phase, 
domain movement does depend on membrane physical 
properties and the radial dependence can be fit by a 
Saffman-Delbriick model with D ~log(l/r). For these 
membranes, we determine 2D viscosities and report acti- 
vation energies for domain diffusion. 
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